Final version: November 10, 2010 



Phobos results on charged particle multiplicity and pseudorapidity distributions in 
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Pseudorapidity distributions of charged particles emitted in Au + Au, Cu + Cu, d + Au, and 
p + p collisions over a wide energy range have been measured using the PHOBOS detector at 
RHIC. The centrality dependence of both the charged particle distributions and the multiplicity at 
midrapidity were measured. Pseudorapidity distributions of charged particles emitted with \r/\ < 5.4, 
which account for between 95% and 99% of the total charged-particle emission associated with 
collision participants, are presented for different collision centralities. Both the midrapidity density, 
dNch/dr], and the total charged-particle multiplicity, N a h, are found to factorize into a product of 
independent functions of collision energy, ^/s NN , and centrality given in terms of the number of 
nucleons participating in the collision, N par t- The total charged particle multiplicity, observed in 
these experiments and those at lower energies, assumes a linear dependence of (lns NJV ) 2 over the 
full range of collision energy of y/s NN =2.7-200 GeV. 

PACS numbers: 25.75.-q, 13.85. Ni, 21.65.+f 



I. INTRODUCTION 

The study of relativistic heavy ion collisions is the only 
known method of creating and studying in the laboratory 
systems with hadronic or partonic degrees of freedom at 
extreme energy and matter density over a significant vol- 
ume. It is for this reason that in recent years such studies 
have attracted much experimental and theoretical inter- 
est, in particular with the likelihood that at the higher 
energies a new state of QCD matter is created. 

During the first five years of the operation of the Rel- 
ativistic Heavy Ion Collider, RHIC, at Brookhaven Na- 
tional Laboratory, the PHOBOS experiment [l[ collected 
extensive data on the production of charged particles over 
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almost the entire solid angle, for a wide range of collision 
energies and colliding nuclei. Many interesting and unex- 
pected results were obtained which have been published 
and their significance discussed in a series of short papers 
una. The early results are summarized and the physics 
interpretation is discussed in Ref. flij . 

This paper presents all PHOBOS results on multi- 
plicity and pseudorapidity distributions, including some 
unpublished data, in a consistent graphical and tabular 
form, together with detailed descriptions of how the re- 
sults were obtained and analyzed. The intention is to 
present the data with a minimum of interpretation. Fit- 
ting of functional forms to the data is done only to fa- 
cilitate reproduction or extrapolation. No significance of 
the functional forms is implied. 

The PHOBOS data cover Au+Au collisions at nucleon- 
nucleon center of mass energy, y/s NN , of 19.6, 56, 62.4, 
130 and 200 GeV, Cu+Cu at 22.4, 62.4 and 200 GeV, 
d+Au at 200 GeV, and p+p at 200 and 410 GeV. Sim- 



ilar measurements, though with less extensive cover- 
age, have been made by the other RHIC experiments 
BRAHMS [13, STAR [16], and PHENIX HJ. These 
measurements extend earlier studies of p + A collisions 
at Fermilab [l8|, , p + A collisions at the Super Proton 
Synchrotron (SPS) at CERN [20], p + Nuclear Emulsion 
[21| . as well as A+A collisions at the SPS reaching ener- 
gies up to ,/s^7 = 17.3 GeV [Hj], and at the Alternating 
Gradient Synchrotron (AGS) at BNL up to 4.9 GeV [H. 
It is expected that heavy-ion collisions will soon be ex- 
tended to higher energies, eventually reaching y/s NN = 
5500 GeV at the Large Hadron Collider at CERN. 

This extensive body of data on the global properties 
of particle production in heavy ion collisions can be used 
to provide insight into both our understanding of the 
mechanisms of particle production and the properties of 
matter that exist at extremes of energy and matter den- 
sities. 

This paper is organized as follows: The PHOBOS 
apparatus relevant for the multiplicity measurements is 
briefly described in Sect. II. This is followed in Sect. Ill 
by a detailed discussion of the data analysis procedure. 
The results are presented in Sect. IV, and a summary is 
given in Sect. V. 



II. EXPERIMENTAL SETUP 

The PHOBOS experiment consists of three major com- 
ponents, a charged particle multiplicity detector covering 
a large fraction of the total solid angle, an array of de- 
tectors for triggering and event characterization, and a 
two-arm magnetic spectrometer used for reconstructing 
the trajectories of a small fraction of the particles emit- 
ted near midrapidity. The entire detector is described in 
greater detail in Ref. Note that only a sub-set of 

detectors were installed for the run resulting in the data 
presented in Ref. [|[ ■ This section will briefly discuss the 
parts of the apparatus used in the current analysis. The 
active areas of several of the detectors used in this work 
as well as the beam pipe are shown in Fig. [T] Note that 
the dimensions of all detectors as well as the positions 
transverse to the beam are shown to scale; the locations 
along the beam have been shifted to facilitate the viewing 
of the detectors in a single figure. The Paddle counter ar- 
ray on one side of the interaction point and the outer four 
Ring counters are not shown. The dimensions and ori- 
entations of the excluded detectors are identical to those 
shown in Fig. [TJ 

The event triggering and centrality determination 
were, for most of the systems, provided by the Paddle de- 
tector, two arrays of 16 plastic scintillator slats positioned 
at ± 3.21 m from the center of the interaction region [24]]. 
Each slat is read out by a single photomultiplier tube con- 
nected to the outer end by a light guide (not shown). The 
active area of the Paddle detectors covers the angular re- 
gion 3 < \r)\ < 4.5, where n = — ln[tan(0/2)] and 9 is 
the polar angle defined with respect to the beam axis z. 
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FIG. 1: The position and orientation of the Be beam pipe and 
the active areas of several of the detectors used in the present 
work. See text for details. 



The primary event trigger required response from at least 
one slat in both counters with a time difference consistent 
with an event occuring near the center of the interaction 
region. Detailed analysis and comparison to simulations 
indicate that this trigger fired for > 97% of the Au+Au 
total nuclear cross-section at y/s~^= 130 and 200 GeV, 
and ~81% for the ^fs^ — 19.6 GeV data [25j], whereas 
the trigger efficiency for Cu + Cu varies between 84% for 
200 GeV and 75% for the 62.4 GeV collisions [13]. The 
same trigger conditions were required for the 200 GeV 
d + Au experiment, resulting in an overall triggering ef- 
ficiency of ~83%, whereas the inelastic p + p collisions 
were obtained by requiring only one slat in one counter 
to trigger in coincidence with the signal from the beam 
bunch crossing clock provided by RHIC [26| . 

The Vertex detector was used in both event characteri- 
zation and multiplicity determination. It consists of four 
layers of Si (silicon) wafers, two above and two below the 
interaction region. This detector covers the two regions, 
each with an azimuthal, 4>, angular extent of roughly 43° 
and r\ range (for events occuring at the center of the in- 
teraction region) of 1 77 1 < 1. The Si detectors are finely 
segmented along the beam direction so that "tracklets" , 
created by combining one hit from the inner and one hit 
from the outer layers, pointed back to the primary inter- 
action point with high accuracy. This vertex location was 
then used to correct the signal in other parts of the mul- 
tiplicity detector (especially the Octagon) for the effect 
of traversing the Si wafers at oblique angles. In addi- 
tion, the distribution of tracklets was used to measure 
the charged particle multiplicity near midrapidity. 

The primary multiplicity detectors are the Octagon 
and the Rings. The former consists of a single layer 
of 92 Si wafers oriented parallel to the beam pipe and 
covering \rj\ < 3.2. Except for regions left open to al- 
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FIG. 2: The geometrical acceptances of the Ring (light), Oc- 
tagon (medium), and Vertex (dark) detectors are shown as 
a function of pseudorapidity, r\, and azimuthal angle, cj>, for 
particles emitted from the nominal interaction point at the 
center of the Octagon array. 



emitted near 77 = 0. As a secondary function, the six 
planes outside the Spectrometer magnetic field, shown in 
Fig. [TJ are also used to measure the charged particle mul- 
tiplicity within the range |Ajj| < 1. Straight-line tracks 
formed by aligning hits from each of the six layers were 
used in this analysis. 

In studying the charged particle multiplicity, it is crit- 
ical that the experimental setup minimizes the modifica- 
tions of the distribution due to absorption, scattering, or 
creation of secondaries. The Be beam pipe is an impor- 
tant element of the PHOBOS design. It is constructed of 
three identical segments, each 4 m long by 76 mm in di- 
ameter, with a wall thickness of m 1 mm. The flanges and 
bolts used to connect the segments are also made of Be. 
The use of a low-Z material and relatively thin walls re- 
duces the interactions of charged particles traversing the 
pipe. The relatively small diameter allows active detector 
elements to be positioned quite close to the interaction 
point as well as at small values of 9 (and hence large rj). 



TABLE I: Geometrical characteristics of Si sensors used in 
charged particle multiplicity measurements. All Si wafers 
have thicknesses in the range 300-340 fim. 



Sensor type 


Active area 


Number of pads 


Pad size 




(mm 2 ) 


Tj x <f> 


(mm 2 ) 


Octagon 


81.28 x 34.88 


30 x4 


2.71x8.71 


Ring 


^3200 


8x8 


« 20-105 


Inner Vtx 


60.58 x 48.18 


256 x 4 


0.47 x 12.04 


Outer Vtx 


60.58 x 48.18 


256 x 2 


0.47 x 24.07 



low unimpeded passage to the Vertex and Spectrometer 
detectors, the Octagon has full azimuthal {(f)) coverage. 
The wafers are segmented in both cf> (about 10°) and rj 
(ranging from 0.06 to 0.005 units depending on distance 
from the center). The Rings consist of six separate de- 
tector arrays (only two are shown in Fig. [T]) located at 
±1.13 m, ±2.35 m, and ±5.05 m along the beam axis, 
extending the coverage for charged particle detection out 
to 1 77 1 < 5.4. These wafers (eight in each Ring) are ori- 
ented perpendicular to the beam and are segmented in 
both <p and 77, with the radial segmentation chosen to 
give approximately constant A77 bin sizes within a single 
detector. These detectors have full <j> coverage. 

The geometrical acceptance of the detectors used in 
the charged particle multiplicity measurements is shown 
in Fig. [21 The pseudorapidity range is here calculated as- 
suming that collisions occur at the center of the octagon 
array. The openings at <f> = 0° and 180° and -0.8< rj <2.1 
are partially covered by the Two-arm spectrometer (not 
shown). Selected characteristics of the Si wafers are listed 
in Table |U The pad dimensions in the longitudinal direc- 
tion are for all Si-sensors sufficient to sort the data into 
77-bins of width A77 < 0.2, which is the bin-size used for 
the multiplicity distributions presented in this work. 

The Si wafers of the Spectrometer detector were de- 
signed primarily to track and identify charged particles 



III. DATA ANALYSIS 
A. Centrality determination 

In interpreting data from heavy ion collisions, the pri- 
mary event characterization parameters are the energy 
of the collision and the overlap of the two nuclei at the 
moment when they interact, commonly referred to as cen- 
trality. Conventionally, ^s NN , the center of mass energy 
available when a single nucleon from one projectile col- 
lides with a single nucleon from the other projectile (ig- 
noring Fermi motion), is used to characterize the energy 
of the collision [3] . 

Characterization of the centrality is more difficult since 
it is not a directly measurable quantity. Conventionally, 
three derived quantities are used as a measure of central- 
ity. The quantity most closely related to measurement 
is the so-called fractional cross-section. For each event, 
the observed energy or multiplicity of charged particles 
emitted into a fixed set of detectors is measured. Events 
are then sorted into bins of fractional total cross-section, 
making the reasonable assumption that the detector re- 
sponse has a monotonic relationship with centrality bins, 
with the bin containing events with the highest multi- 
plicity corresponding to the most central collisions, i.e. 
those with the smallest impact parameter. 

In PHOBOS, for the Au+Au data obtained at y/s NN 
= 56, 62.4, 130, and 200 GeV, and the Cu+Cu data 
at 62.4 and 200 GeV, it is the energy deposited in the 
Paddle detector arrays [24J that is used for sorting the 
data into bins of fractional cross-section. Since the trig- 
gering efficiency plays an important role, especially for 
low-multiplicity peripheral events, it is important to ac- 
curately estimate this quantity. This is done by de- 
tailed comparisons with Monte Carlo simulations as de- 
scribed in detail in Ref. [25| . The Paddle detectors cover 
the range 3< \r]\ <4.5, into which a large multiplicity 
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of charged particles are emitted at these energies. For 
ys^7=19.6 GeV Au+Au and 22.4 GeV Cu+Cu colli- 
sions, however, the charged particle multiplicity is small 
in the large \r\\ region covered by the trigger detectors, 
such that these centrality measurements become less re- 
liable. At these lower energies a different measure, the 
total energy deposited in the Octagon detector in the 
range -3.0< r\ <3.0, is used for sorting the data. The 
centrality measurements for d+Au collisions represent a 
special challenge. For reasons discussed in Ref.[l(3|, the 
six Ring detectors spanning 3.0< \r)\ <5.4 were used to 
obtain the most consistent centrality measurement. 

A detailed discussion of the choices of these detectors 
and of the fractional cross-section determination in PHO- 
BOS can be found in Ref. [lj]. 

There are two other measures of centrality that have 
been found more useful when it comes to comparison of 
heavy- ion data with p+p data, as well as for the inter- 
pretation of the data, and which can be derived from the 
fractional cross-section. Both are motivated by the fact 
that, because of relativistic time dilation, the collision 
duration time at high energies is very short compared to 
the typical time-scale for soft particle production and for 
nuclear rearrangement or movement of nucleons within 
the nuclei. Assuming that the nucleons in each nucleus 
are indeed unaltered and fixed in the transverse direction 
during the time of the collision of the two nuclei, one can 
introduce the concept of Neon, the number of individual 
nucleon-nucleon collisions that occur during the nucleus- 
nucleus collision. One can also introduce the concept of 
Np ar t, the number of nucleons that have made at least 
one collision with a nucleon of the other nucleus dur- 
ing the collision. The latter is the same quantity as the 
number of "wounded nucleons" introduced by Bialas et 
al. [27[ to interpret the p+A results obtained by Fermilab 
experiment El 78 [28} . 

Both Np ar t and N co u can be derived from the frac- 
tional cross-section by modeling the collisions of the two 
nuclei and assuming that the fractional cross-section in- 
creases monotonically with N part or N co u. In PHOBOS, 
we model the collisions assuming the nuclei are a collec- 
tion of hard spheres, with radii corresponding to the in- 
elastic nucleon-nucleon cross-section at the appropriate 
energy, distributed according to the Wood-Saxon func- 
tional form for Cu and Au ions and the Hulthen wave 
function for deuterons, and that the nucleons are unal- 
tered and follow straight-line trajectories as they make 
collisions inside the nucleus (a procedure often referred 
to as the Glauber model (29[). Details of the modeling 
can be found in Ref 

Figure [3] shows the total energy deposition spectrum 
from the Paddle detectors covering 3< |?7| <4.5 as well as 
the deduced number of participants, N part , for the three 
centrality bins of 0-3%, 20-25%, and 45-50% of the total 
cross section in 200 GeV Au + Au collisions (indicated 
by shaded bands in Fig. [3^,). The corresponding N part 
distributions are shown in Fig. \3jp and values of the cen- 
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FIG. 3: Illustration of the conversion from paddle signal to 
Np a rt for 200 GeV Au + Au collisions. Panel a) shows the 
ranges of the Paddle signals corresponding to the 0-3%, 20- 
25%, and 45-50% most central collision. The corresponding 
distributions of the number of participants, N part , are shown 
in panel (b). 



cross section bins are listed in Table [V] (see Appendix) . 



B. Signal processing 

Nine types of silicon pad sensors are used in PHO- 
BOS. These sensors are segmented into as many as 1536 
individual pads in the case of first spectrometer detector 
plane. Each pad is connected to one channel of a 128- 
or 64-channel pre-amplifier and readout chip. When a 
particle traverses the detector a trigger-derived "hold" 
signal causes a front-end chip to store the signal peak on 
all channels simultaneously and later multiplexes them 
onto a differential analog output bus. The technical de- 
tails of the Si-detectors and the associated electronics are 
described in Ref. [l[ and references therein [30l - l33l |. 

For every event satisfying the trigger criteria, the sil- 
icon data are converted from raw 12-bit ADC values to 
calibrated deposited energies by applying a series of four 
algorithms, namely: 1) pedestal subtraction, 2) com- 
mon mode noise correction, 3) gain correction, and 4) 
hit merging. 

The pedestals of the read-out system are subtracted 
from the raw signal data event-by-event for each channel. 
The pedestal widths provide a measure of the noise in the 
system that arises due to silicon detector leakage cur- 
rents, as well as intrinsic electronics noise. The pedestals 
are obtained from dedicated runs with triggers derived 
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FIG. 4: a) Energy deposition spectra for 2.5 < r\ < 3.2 for the 
Octagon detector, b) Same as a) but for 77 < 1. The dotted 
histogram corresponds to the raw spectrum, whereas the solid 
histogram shows the effect of the signal merging correction. 




u 

FIG. 5: (color online) (a) Map of the raw energy deposition 
versus pseudorapidity r\ for the Octagon detector, (b) Same as 
(a) after applying the merging of signals from particle travers- 
ing multiple pads and correcting for the angle of incidence. 



from a random electronic pulser with adjustments using 
low multiplicity collision data. 

The detector readout system also has a noise compo- 
nent that fluctuates with time, but is common to several 
channels. The level of this "common mode" noise (CMN) 
is estimated on an event- by- event basis. The most prob- 
able value within a certain range of the pedestal sub- 
tracted ADC signals of all channels from one read-out 
chip is considered to represent the CMN for the chip. 
This value is subtracted from all channels of the relevant 
chip. 

For very high occupancies, the data-based common 
mode noise correction in the octagon detector (|?7| < 
3.2) becomes slightly inaccurate since it relies on hav- 
ing enough empty channels to provide a baseline. The 
more highly segmented vertex detector allows us to mea- 
sure this effect and make a further data-based correction. 
This further correction was only required near midra- 
pidity (\eta\ < 1.5) for the central (0-10%), high energy 
(y/s NN > ISOGeV) Au+Au data, and it was less than 
4% everywhere. 

The conversion from ADC channel to energy deposi- 
tion in the Si-wafer is obtained by using a calibration 
circuit 33], which injects a known amount of charge into 
the front end of each channel during special calibration 
runs. We have observed that the gain of the front-end 
chips is linear over most of the range. Non-linearities are 
observed only for very large signals corresponding to the 



energy deposition of around 65 minimum ionizing par- 
ticles (65 MIPs), where 1 MIP corresponds to about 80 
keV most probable energy loss in 300 silicon. 

In the center of the Octagon and in the Rings, pri- 
mary particles were incident very close to normal to the 
silicon wafers (depending on the location of the primary 
interaction), and, therefore, deposited an amount of en- 
ergy corresponding to that of a minimum ionizing parti- 
cle. Near the Octagon extremities, a particle's trajectory 
is at very oblique angles, and the particle may deposit 
energy in two or more adjacent pads. In this case the 
deposited energy from all pads contributing to a single 
particle track was merged. The location of each pad with 
a deposited energy above the hit threshold of 
keV was identified. For pads where the detector geome- 
try indicates that the particle track must traverse more 
than a single pad, the nearest neighbor pads are checked 
to determine if they also contain some deposited energy. 
These signals are added together, and assigned a pseudo- 
rapidity value appropriate for the pad that contained the 
largest energy signal for that track. The nearest-neighbor 
pads are no longer considered in the subsequent analysis. 
FigureH^, illustrates the effects of this merging procedure. 
The dashed and solid histograms show deposited-energy 
spectra for unmerged and merged hits, respectively, over 
an 77 range of 2.5 to 3.2. The improvement in the signal 
quality is evident. This correction is much smaller near 
midrapidity as illustrated in Fig. HJd, which corresponds 
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to 1 77 1 < 1. Note that these spectra have already been 
corrected for the angle-of-incidence as described below. 

After merging hits, the deposited energy, AE tr , was 
corrected for its angle of incidence, as well as the thick- 
ness of the silicon wafer, so that the deposited energy 
could be compared to a common most probable value 
of AEmip=82 keV. The uncorrected energy deposition 
(normalized to a common wafer thickness of 300 /im; the 
thickness of the Si wafers were measured and found to 
lie in the range 300-340 /jm) is shown as a function of 
r] in FigJS^,. For valid primary tracks, we observe that 
the energy deposition for a single minimum ionizing par- 
ticle follows a AE fa cosh r\ dependence corresponding to 
the length of the track in the Si wafer. Figure [Sja shows 
the angle-of-incidence corrected energy deposition, AE, 
plotted as a function of 77. To select primary tracks cre- 
ated at the collision vertex, only tracks with AE values 
greater than 0.6xA_Em/p= 45 keV were included in the 
subsequent analysis. This selection eliminated a large 
fraction of tracks from secondary sources, e.g. the beam- 
pipe, magnet yokes, etc., that have a small AE value 
inconsistent with primary particles. Many such tracks 
are seen in the region of small AE values and large \r)\ 
in FigfSjD. This background was insignificant for small 
values of r\, where the amount of extra material between 
the interaction point and the detectors is minimal. 

Malfunctioning channels are flagged and stored in a 
"dead channel" map. Additional dead channels are found 
by looping over events in a reference run. When the 
signal in a channel was 5 times larger than its noise it was 
considered a valid hit and a hit-counter was increased 
and the energy value added to an energy sum. At the 
end of a run, the average energy a particle deposited in 
a silicon pad was calculated and compared to a set of 
criteria, which included a minimum and maximum hit 
occupancy for the channel and a range of the average 
energy deposition per hit. Channels that failed to satisfy 
these criteria were marked "dead" in a map. 



C. Vertex finding 

At RHIC, the ion beams collide at zero degrees, such 
that the interaction vertices are distributed over a rel- 
atively large region of about 1 m along the beam-line 
around the nominal collision point. An approximate 
measurement (±7.5 cm) of the z-coordinate of the colli- 
sion vertex was obtained from the relative time difference 
between the paddle counters used in the event triggering. 
Also, the density distribution along the z-axis of hits with 
an energy deposition between 0.5 and 2 AEmip yielded 
a vertex finding accuracy of Az ~ 0.5 cm. Although 
these detectors provide relatively poor vertex determina- 
tion, they are effective for peripheral collisions and give 
an important cross check on the vertex reconstruction 
obtained from the Vertex detector and Spectrometers. 

The Vertex detector was effective in determining the 
z and y coordinates of the collision vertex by connecting 



hits in the outer and inner silicon layers to form "track- 
lets" - two point tracks. All possible combinations of hits 
on the inner and outer layer of the upper and lower Ver- 
tex detectors are constructed, producing a distribution 
of all possible vertex positions. When the tracklets point 
to the correct vertex position, a peak is observed in the 
vertex position distribution, whereas incorrect combina- 
tions contribute to a broad combinatorial background. 
The Vertex detector can only determine the x coordi- 
nate with an error much larger than the beam profile 
because of the large dimension (24 mm and 12 mm) of 
the pads in that direction. Instead, this x coordinate 
was accurately determined by using straight-line tracks 
formed on the basis of hits in the first 6 silicon layers of 
the Spectrometer arms. 

The final vertex location was determined by an algo- 
rithm that performed an arbitration between the various 
measurements. An overall accuracy of a x ~ 0.15 mm, 
Oy ~ 0.15 mm, and a z ~ 0.06 mm was obtained for the 
15% most central Au + Au collisions at v /s^7=200 GeV 
[34| . Somewhat poorer resolutions were achieved for the 
lower multiplicity peripheral events and at lower collision 
energies. The efficiency of the vertex reconstruction has 
been evaluated on the basis of Monte Carlo simulations 
using events with vertices in the range \z\ <10 cm. In the 
case of the vertex detector method it was greater than 
85% for the 40-45% centrality bin increasing to 100% for 
the 30% most central collisions at y^s NN —130 and 200 
GeV. Because of the lower multiplicity at - v /s JVJV =19.6 
GeV, the vertex reconstruction efficiency was reduced in 
this case. 

For the d+Au and p+p collisions, the small number 
of tracks in the Vertex detector renders the precise Ver- 
tex detector and Spectrometer methods too inefficient for 
practical use. A new method, based on the energy depo- 
sition in adjacent detector pads in the Octagon detector, 
retains high efficiency even for low multiplicity events, 
but it results in a less accurate determination of the ver- 
tex position in the range of a z — 0.5 to 2.0 cm depending 
on the event multiplicity, see Ref. [35[ for details. Even 
this approximate determination of the vertex position is 
sufficient for the extraction of the charged-particle multi- 
plicity. The reconstructed multiplicity of an event is only 
slightly affected by the vertex position error while the 
more significant modifications of the shape of dN^/dr] 
mostly cancel by averaging over many events. 



D. Monte Carlo simulations 

Although the reconstruction of the event multiplicity 
eliminates most of the background on the basis of energy 
deposition in the high-77 region of the Octagon detector 
and by accepting only tracks that point back to the ver- 
tex position in the two-layer Vertex detector, some back- 
ground hits remain. This residual background is esti- 
mated on the basis of Monte Carlo simulations and ap- 
propriate corrections are then applied. Such simulations 



0.2- . 
- Ring 

0.0 L 



Ring 



-200 



-100 



Octagon 

_^ i , t 



z (cm) 



Ring 
100 



Ring 



200 



FIG. 6: a) Map of simulated secondary particle sources near 
the horizontal plane, y=0. For clarity, interactions in the air 
were neglected. The beam pipe with flanges, silicon sensors 
and elements of their support structures are visible, b) The 
ratio of hits from secondary particles to all hits in the multi- 
plicity sensors located at different positions along the z-axis. 
The values are averaged over the azimuthal angle <j>. The 
farthest Ring counters are not shown. 



FIG. 7: Illustration of tracklet counting using two consecutive 
Si layers, in this case Vertex detectors. Hits in the inner and 
outer Si layers are connected to form "tracklets" that point 
back to the vertex position in the center of the figure. 



at z=±5.05 m are identical at 0.545. 



of the detector response and sources of background parti- 
cles, that do not originate from the collision vertex, were 
carried out using the HIJING (3(| event generator. The 
GEANT 3.21 [37[ package was used to simulate the de- 
tector response and particle tracking using the detailed 
geometry of the PHOBOS detector. Simulations using 
events generated with the VENUS [H and RQMD 
codes were also performed to estimate systematic errors. 
All elements of the sensors and support structures were 
precisely modeled, especially those near the interaction 
point and along the path of particles to the more distant 
sensitive elements. In order to illustrate the degree of 
details of the background simulation in the PHOBOS ap- 
paratus we present in Fig. [5] a map of the origin points of 
secondary particles leaving hits in the multiplicity detec- 
tors. One observes that even small elements, such as the 
Si sensors of the Spectrometer, appear in this map. Note 
that this figure only includes secondary particles gener- 
ated near the horizontal plane. Therefore it does not 
show particles created in much more massive elements, 
such as the magnet yokes, which are located above and 
below the horizontal plane. 

Using Monte Carlo studies, the fraction of hits from 
secondary particles to the total number of hits in the 
multiplicity detectors was determined; this is shown in 
Fig. [U>. This fraction was smallest for hits in the Oc- 
tagon sensors closest to the interaction point but it in- 
creases towards the extremities of the Octagon array be- 
cause of the smaller solid angle with respect to the vertex 
position subtended per unit area. The ratio decreases for 
the Ring counters. The asymmetry with respect to r]=0 
was caused by secondary particles created by the mag- 
net yokes, located in the positive z region. The values of 
background fraction for the two Ring detectors located 



E. Tracklet counting 

Three different methods were employed to determine 
the charged particle multiplicity. The "tracklet count- 
ing" method was used in the (rj, </>) range covered by two 
or more layers of silicon detectors, namely the regions 
covered by the Vertex detector and the Spectrometers, 
where it is possible to correlate "hits" in the two lay- 
ers and thereby construct "tracklets" that point back 
to a previously determined vertex position. An illus- 
tration of tracklets originating from the vertex position 
is shown in Fig. [71 The multiplicity measurements pre- 
sented here were performed without the PHOBOS mag- 
net being powered so that straight-line tracks may be 
assumed except for the effects of small angle scattering 
in the Be beam pipe and the intervening silicon layers. 

The Vertex detector was used in the tracklet analysis 
for particles emitted in range \r]\ <1. The tracklet recon- 
struction procedure generated "seed" tracks using hits in 
one silicon layer and the reconstructed vertex position. 
A seed track was extrapolated to the other silicon layer 
(the "search" layer) and the location where the seed track 
traversed the "search" layer was compared with "hits" in 
that layer. The conditions 

\5(f>\ = \(f) see d - 4>search\ < 0.3 (1) 
\Sf]\ = \r/seed ~ Vsearch\ < 0.1 

were required for a tracklet candidate ■ This method 
allowed a single hit in the "search" layer to be shared 
between more than one tracklet. In such cases only one 
tracklet was accepted for further analysis. The combi- 
natorial background was estimated by performing the 
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FIG. 8: Illustration of the combinatorial background sub- 
traction for non-vertex tracks for a) Vertex detector and b) 
Spectrometer tracklet analysis. The solid points represent 
the actual detector geometry, whereas the open points were 
obtained after a 180° rotation of one of the detector layers 
around the beam axis in order to estimate the combinatorial 
background in the acceptance region represented by vertical 
dashed-dotted lines. 



same tracklet-finding procedure with the "search" layer 
rotated by 180° in <fi, see Fig. [5^,. After subtraction of 
the combinatorial background, the charged-particle mul- 
tiplicity was computed taking into account the centrality 
and z vtx dependent detection efficiencies, which were ob- 
tained from Monte Carlo simulations. 

The tracklet analysis method was also applied to the 
first and second layers as well as the fifth and sixth layers 
of the silicon detectors in the two-arm Spectrometer, see 
Fig-Eb- This allowed the multiplicity to be determined in 
the < rj < 1 range. In this method, a slightly different 
tracklet acceptance criterion was used, namely 



V^^W < 0.016 



(2) 



instead of the "rectangular" acceptance region described 
by Eq. [TJ Tracklet multiplicities obtained from the two 
detectors were in excellent agreement (typically <5%). 

Tracklets were used to obtain the centrality depen- 
dence of the charged particle multiplicity in the midra- 
pidity region, \r]\ < 1, and to validate measurements ob- 
tained from the hit counting and energy deposition meth- 
ods described in the following section, which were applied 
to the full ry-range subtended by single-layer silicon de- 
tectors (Octagon and Ring detectors). 



F. Single Si layer analysis 

The energy deposited in the single layer of silicon de- 
tectors making up the Octagon and Ring multiplicity de- 
tectors was used to determine the charged-particle mul- 
tiplicity over the entire r] range accessible to PHOBOS. 
The energy deposited in these detectors for a typical cen- 
tral event is shown in FiglHl 

The charged-particle multiplicity density was deter- 
mined from the deposited energy in two independent 
ways. In both cases, only events with a reconstructed 
vertex on the beam axis within ±10 cm of the geometri- 
cal center of the multiplicity detector were used. The two 
methods were somewhat complementary. Hit counting 
was relatively insensitive to the energy calibration since 
it relied only on whether the signal in a silicon pad was 
above a certain detection threshold. On the other hand, 
the energy deposition method used more of the avail- 
able information at the cost of increased dependence on 
Monte Carlo simulations by associating the signal height 
in a pad with an estimated number of traversing charged 
particles. This method is therefore directly sensitive to 
the energy calibration. Both methods are described in 
the following sub-sections. 



G. Hit counting method 

This method was based on counting pads with an en- 
ergy deposition greater than a pre-determined minimum 
value. The charged-particle multiplicity was determined 
by summing over all pads lying in a pseudorapidity range, 
77, for events within a certain centrality, b. The pseudora- 
pidity density dN c h / drj was obtained using the expression 



1 



dN r 



drj N(b) 



N(b) 

E 



JVjwt. 



O(M) x /fefcg(M) 
A{z,rj) x Arj 



(3) 



Here, N(b) is the total number of events for a given 
centrality b, Nhits is the number of hit pads, 0(b, if) rep- 
resents a correction due to detector occupancy, A(z,rf) 
is the geometrical detector acceptance, which varies with 
both the position of the collision vertex, z, and the pseu- 
dorapidity, and the quantity ftkg takes into account the 
effects of backgrounds and secondary particle production. 
The determination of these correction terms is described 
below. 



1. Occupancy 

Having corrected the deposited track energies for their 
angle of incidence, and eliminated low-AE? tracks orig- 
inating from background and secondary particles, the 
data must be corrected for the effects of detector occu- 
pancy. In particular, for central events over the entire 
rj range and for non-central events in the \rf\ < 3 range 
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FIG. 9: (Color online) Single event energy deposition map in 
cj> vs. rj for a central Au+Au collision at 130 GeV. The open 
areas at cj>=0 and <f> — n are openings in the Octagon detector 
in front of the Spectrometers. Small open areas at <j> — ±7r/2 
and r] about -2 and 1 in front of the Vertex detectors are also 
visible. 
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there is a significant probability that more than one par- 
ticle had traversed a single detector pad. The detector 
occupancy, 0(6, 77), depends upon both 6 and 77 and was 
determined using two different methods. 

The first method relies on the assumption that, for a 
given centrality, 6, and pseudorapidity range, A77, the 
probability that N particles pass through a pad is de- 
scribed by a Poisson-statistical distribution, 



P(N) = 



N\ 



(4) 



where (i is the average number of tracks per pad deter- 
mined over an ensemble of events. The value of \x depends 
on both centrality and pseudorapidity, and was deter- 
mined experimentally from the event sample. To deter- 
mine 6), the number of pads, Nhu, in a A77 range 
with a valid energy signal was compared to the number 
of pads with no hits, N nohit . The ratio R = N hlt /N nohlt 
is related to /j, by n(rj,b) = ln[l + R(r],b)}. The occu- 
pancy correction factor is then given by Poisson statistics 
0(r],b) = n(rj,b)/ (exp fJ.(r],b) — 1). Typical values for the 
occupancy correction factor were between 1.0 and 1.2 for 
large values of 77, or for peripheral collisions, but were 
larger (0(77,6) ps 1.8) for central Au + Au collisions at 
i] s» at a collision energy of ^s NN — 200 GeV, see 
Fig. HJfti. 

An independent experimental confirmation of the va- 
lidity of the Poisson occupancy determination was made 
by studies of the AE pac [ spectra as a function of 77 and 
6. After merging hits, and correcting the resulting track 
energies for angle of incidence, the AE pa d spectra showed 
structures characteristic of the energy deposited by one, 
two, or more particles traversing a single pad (see Figj4j). 
Each AEp a( i spectrum was fit to a set of Landau func- 
tions convoluted with Gaussian functions to account for 
the intrinsic energy resolution of the pad. From the re- 



FIG. 10: (a) The single pad occupancy correction factor 0(r/) 
for central (0-6%) and peripheral (45-55%) Au+Au collisions 
at 130 GeV are shown as solid and open symbols, respec- 
tively. The data for the Octagon detector are shown as circles, 
whereas the triangles, diamonds, and squares correspond to 
Ring detectors at increasing distances from the vertex posi- 
tion. The background correction fbkg{rj) factor and the total 
correction factor 0(77) x fbk g {rj) are shown in panels (b) and 
(c), respectively, using the same symbols. 



suits of these fits, the relative contributions of N =1, 2, 3 
hits to the AE pa d spectra could be determined as a func- 
tion of 6 and 77. These relative contributions then yielded 
an independent occupancy correction factor 0(6, 77). This 
method of determining the detector occupancy was used 
as a confirmation of the validity of the Poisson occupancy 
method. 



2. Background and Monte-Carlo Corrections 

The majority of particles produced by secondary or 
background interactions were eliminated by requiring the 
value of AE pa d to exceed a threshold value as described 
above. However, there still exist additional background 
contributions that could not be eliminated by using mea- 
sured quantities alone. Such backgrounds included sec- 
ondary particles produced in the beam pipe, in the mag- 
net yoke, or in other detector elements in addition to 
those generated via feed-down from the weak decays. The 
background correction also accounts for the absorption of 
particles in the beam pipe before they had reached the 
silicon detectors. To account for these effects, Monte 
Carlo simulations were used to determine the response of 
the apparatus to particles produced by a variety of event 
generators, as described in Sect. IIIID1 Such simulated 
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event data, generated for a variety of collision central- 
ities were analyzed in exactly the same manner as real 
data to estimate the residual background effects. 

The simulated data were first passed through the al- 
gorithm used to determine 0(b,r)) following the same 
method as for the experimental data. The simulated data 
were subsequently corrected for detector acceptance ac- 
cording to Eqf51 without the term fbkg- The resulting 
dNch/dr] generally did not agree with the "true" known 
form of dN c h/drj obtained directly from the output of the 
event generator. The final correction factors that were 
applied to the experimental data were determined from 
a comparison of the reconstructed Monte-Carlo dN c h/drj 
to the "true" distributions as 



hkg(b,ri) 



dN ch dN ch 

-(true)/ — — (recon) 



drj 



dq 



(5) 



Examples of these correction factors are shown as a 
function of 77, for both peripheral, and central collisions, 
in Fig. [TDb . As expected, near 77 = 0, fbkg ~ 1 indi- 
cating that there is only a small contribution from these 
processes. Near r/=+3 to 3.5, however, fbkg is somewhat 
smaller due to a large number of secondary particles pro- 
duced in the steel of the magnet yoke. As expected, the 
number of background and secondary particles, as well 
as the number of particles absorbed in the beam pipe is 
proportional to the number of primary particles. This 
can be seen from the fact that fbkg is almost identical for 
peripheral and central collisions (see Fig. HOb ) . 



H. Energy deposition method 

A second, largely independent, method of determining 
dN c h/drj used the AE pac i signal to estimate the number 
of tracks traversing a pad. This method uses the energy 
deposition in the first layer of all the silicon detectors 
in PHOBOS, including not only the Octagon and Ring 
multiplicity detectors, but also the inner layers of the 
Vertex detectors, and the first Spectrometer layers. The 
"analog" multiplicity method obtains the total dN c h/drj 
from the relation 



dNci 
d-q 




AEp ac [ X fp r i X f a b 

Ar] pad x AEtrk 



(6) 



where the sum extends over all active pads in the detec- 
tors used for the analog analysis. AE pac i is the energy 
deposited in a single pad, AE tr k is the average energy 
deposited in 300 /im of silicon by a single ionizing parti- 
cle (SIP), fpriiv) is the fraction of primary particle tracks 
out of the total number of tracks in a given pad, f a b s is a 
correction that takes into account the effects of absorp- 
tion in the beam pipe, and Ar] pa d is the acceptance in 77 
for a given pad. The quantities AE tr k and f pr i depend 
on the relative position of a given detector sensor and 
the collision vertex as described below. The deposited 
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FIG. 11: (a) The energy deposition per detector pad AE pa d, 
(b) the energy deposition per particle track AE tr k, (c) the 
fraction of primary tracks f pr j , and (d) the absorption cor- 
rection factor f a i, s are shown as a function of r\ for central 
(0-6%) (solid symbols) and peripheral Au+Au collisions (45- 
55%) (open symbols) at 130 GeV. The circles and squares 
refer to simulations for the Octagon and Ring detectors, re- 
spectively 



energy per pad AE pa d is corrected for the angle of in- 
cidence based on the measured vertex position. Figure 
1111 shows the pseudorapidity dependence of the various 
quantities in Eq. 6 for 200 GeV central (0-6%) and pe- 
ripheral (45-55%) Au + Au collisions, and that the ratio 
f P ri/ AEtrk used in the calculation of dN c h/dr] is nearly 
independent of rj. 



1. Analog method correction parameters 

The average deposited energy per track, AE tr k, is 
determined from Geant simulations (as described in 
Sect. IIII Dp . The value of AE tr k ranges between 50 
and 160 keV, depending on the direction of the track 
traversing the sensor and the mixture of primary and 
secondary particles. Note that the average energy depo- 
sition, AE t rk , is substantially larger than the most prob- 
able value clue to the asymmetric nature of the Landau 
distribution as illustrated in Fig. 2] The r\ dependence in 
the Octagon primarily reflects a sensitivity to the broad- 
ening of the energy-loss distribution as the particles en- 
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counter more silicon with increasingly shallow angles of 
incidence. The additional ^-dependence is due to the 
changing mix of particle types, particle momenta, and 
the number of secondary particles produced within the 
volume of the silicon. These variations do not depend 
strongly on the event generator used to obtain the pri- 
mary particle distributions. The values obtained here are 
to be compared with the mean energy deposition value 
of 117 keV noted in the Particle Data Book[41] for min- 
imum ionizing particles normally incident on a 300 /im 
thick wafer of silicon. The value of AE tr j. is calculated 
for each sensor in the first silicon layer in the PHOBOS 
detector, and its dependence on r\ is fitted with a poly- 
nomial parameterization for computational purposes. 

The quantity, AE tr k, was also measured using the ac- 
tual data under low occupancy conditions, corrected for 
multiple occupancy. The final validation of the analog 
method parameters, and the systematic uncertainty as- 
signments, are provided by observing the close agreement 
between the two (and three, where available) analysis 
methods. 

The fraction of primary particles to the total number 
of particles f pr i = N ^l^i^ is determined from simula- 
tions using events produced by HIJING (see Sect. IHI 01 . 
and includes secondary particles produced by interac- 
tions with material as well as feed-down from weak de- 
cays. The values of this parameter vary from approx- 
imately 0.9 near 77=0, to ~ 0.3 near the ends of the 
Octagon, or for the Ring detectors closest to the vertex 
position. The dependence of f pr i on pseudorapidity, cen- 
trality and vertex position is similar to that of the back- 
ground correction factors from the hit counting analysis 
described above. This quantity is also parameterized for 
computational purposes. 

The absorption coefficient f a b s takes into account the 
absorption of particles in the material traversed before 
encountering the silicon detectors, chiefly the 1 mm Be 
beam pipe. As such, this correction has a cosh?y depen- 
dence, with the value of f a bs being approximately 0.98 
near 77=0. The acceptance of each pad, A77, is calculated 
as a function of rj from the measured PHOBOS geometry. 
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FIG. 12: The energy dependence of the participant-scaled 
charged particle multiplicity dN c h/dri/ '(N par t/2) at midrapid- 
ity \v\ < 1 is shown for 0-6% central nucleus-nucleus collisions 
(solid symbols) and compared to pp/pp non-single diffractive 
(NSD) (open diamonds) [43 - 4451 ] and inelastic (open squares) 
[424441 . 14a44a |) and (crosses) present work. The solid squares 
and triangles represent results from the PHOBOS experiment 
for Au + Au and Cu + Cu collisions, respectively. The solid 
circles and solid diamonds are obtained from SPS Pb + Pb 
[22] ] and AGS Au + Au Hi] collisions, respectively. The 56 
GeV Au + Au data point is from Ref. Q]. The solid line is 
a linear fit to the Au + Au and Pb + Pb data, whereas the 
dashed curve is a fit to the inelastic pp/pp data. See text for 
details. 



A. Charged particle production at midrapidity. 

1. Energy dependence 



IV. RESULTS 

The results are for all charged particles excluding weak 
decays and corrected for missing low px particles. Al- 
though the tracklet and single-layer measurements of the 
charged particle multiplicity are in good agreement, we 
present the results separately, because of their different 
range of applicability. For the centrality and energy de- 
pendence of the charged particle density at midrapidity 
the tracklet method is more accurate, whereas the over- 
all pseudorapidity distribution can be obtained only from 
the hit-counting and energy deposition analysis which 
can be applied away from 77 = 0. 



The energy dependence of the charged particle multi- 
plicity, normalized by the number of participant pairs, 
is shown in Fig. [T2b for central (0-6%) Au + Au (solid 
squares and diamonds), Cu+Cu (solid triangles), Pb+Pb 
(solid circles) [22j collisions. The Au + Au and Cu + Cu 
data shown in this figure and listed in Table [TT] rep- 
resent an average of the results from the tracklet and 
single Si layer analyzes. The heavy-ion data are com- 
pared to the dN c h/drj data obtained in non-single diffrac- 
tive (open diamonds) [HI, [H, [45[, and inelastic (open 
squares [H, [H, Hgj]) p + p or p + p-collisions. One ob- 
serves a strong enhancement of the midrapidity charged 
particle production per participant pair in Au + Au col- 
lisions compared to p + p of up to 70%. 

The overall energy dependence of jf^—^^f is loga- 
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FIG. 13: Panel a: The charged particle multiplicity at midra- 
pidity \n\ < 1 per participant pair, dN c h/dr}/ '{N par t /2) ob- 
tained from the "tracklet" analysis is shown as a function of 
{Npart) (solid symbols) for Au + Au collisions. The shaded 
ovals indicate the 90% confidence limit systematic errors. For 
comparison, open points show the "single Si layer" analysis. 
The solid curves represent fits to the data using the form 
given in Eq. [8]- [10] (excluding the pp points). The pp points 
at (N par t)=2 were interpolated using the fit to the inelastic 
pp/pp data, see Fig. 1121 See text for a discussion of panel (b). 



rithmic for A+A collisions. The solid line, given by 
2 dN ch 



part/ 



(hj 



= 0.781n(Vs) - 0.4, 



(7) 



is seen to describe the Au + Au and Pb + Pb data quite 
accurately over the two orders of magnitude of collision 
energy, but it also appears that the Cu + Cu data fall 
slightly below this trend. Also the 56 and 62.4 GeV 
Au + Au points fall slightly below this line, which may 
indicate a curvature to the collision energy dependence. 
The dashed curve is a fit to the inelastic pp and pp data, 
namely dN ch /dr]/{N par t/2) = 0.35 + 0.52s - 123 . 



2. Centrality dependence and factorization 

The midrapidity charged-particle multiplicities nor- 
malized to the number of participant pairs, (N part ) /2, are 
shown for Au + Au collisions (solid symbols) in Fig. [FJ] 
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FIG. 14: Same as Fig. 1131 but for Cu + Cu collisions. 



TABLE II: Summary of the midrapidity ||,,|<i/(A r P art/2} 
charged particle multiplicity for Au+Au and Cu+Cu for 0-6% 
central collisions. The data are averaged over those obtained 
from the tracklet counting and single Si layer analysis. Errors 
represent averages of the 90% C.L. systematic errors for the 
two methods. Statistical errors are negligible. 



System , 


/S^(GeV) {Npart} 


dN ch/ dr l\\r ] \ < l 

{N mTt /2) 


Au+Au 


19.6 


337 


1.87+0.15 




56 


330 


2.47+0.27 




62.4 


338 


2.64+0.20 


» 


130 


342 


3.35+0.25 




200 


345 


3.76+0.33 


Cu+Cu 

55 


22.4 
62.4 
200 


99 
96 
100 


1.94+0.15 
2.47+0.19 
3.48+0.28 



as a function of centrality of the collision expressed by 
{Np ar t) and listed in Table fVl column 4 (see Appendix). 
The corresponding data for Cu + Cu collisions are given 
in Fig. [Hi and Table ED (see Appendix). The shaded 
ovals represent estimates of the 90% confidence level sys- 
tematic errors in the measured dN c h/dr) values and the 
calculated number of participant pairs, (N part )/2. One 
observes that the normalized charged particle produc- 
tion at all energies increases with (N part ) and exceeds 



13 



the measurements for pp and pp inelastic collisions repre- 
sented by stars and open triangles at N part =2 in Figs. [TBI 
and [U The 200 GeV pp value of 2.25±0.1 was mea- 
sured in Ref (43[, whereas in 62.4 GeV collisions a value of 
1.89±0.1 was measured in Ref. [42j]. A fit to the multiplic- 
ities observed for inelastic pp collisions Ref. [13, EJ S| 
was obtained, see Fig. rpJl and used to derive data points 
for 19.6, 22.4, and 130 GeV collisions. 

As demonstrated in a previous PHOBOS publica- 
tion |9j, the collision energy and centrality dependences 
of charged particle production in Au + Au collision at 
midrapidity exhibit factorization such that 



2 dN c h\\ v \<i 



part/ 



d/q 



= f(s) x g(N part ). 



(8) 



For Au + Au collisions we find that the data are very 
well described by the functions 



f(s) = 0.0147(lns) 2 +0.6 
g(N part ) = l + 0.095iV : 



1/3 
part 



(9) 
(10) 



shown as solid curves in Fig. [T3k . 

The ratio of the data to this fit is shown in Fig. [J3b. 
The small standard deviation, a — 0.0155 (shaded band) 
from the mean value (data/ fit) =0.9993 of all Au + Au 
data points (horizontal line) illustrates the accuracy of 
the factorization; a is much smaller than the estimated 
total error (systematic and statistical) on the data points. 
Note that the functional form of the energy dependence 
chosen here is different from the overall trend discussed 
in the previous sub-section since it provides a slightly 
better fit over the limited range of collision energy for 
the Au + Au data shown in Fig. [T2J this choice illustrates 
better the high degree of energy-centrality factorization 
observed in these data. 

For Cu + Cu collisions the same energy dependence 
function /(s) applies, whereas the N part dependence is 
given by 
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FIG. 15: dN c h/dri vs r\ (solid points) for ten centrality bins 
representing 45% of the total cross section for Au + Au col- 
lisions at ^/s NN =19.6 GeV. The solid curve is a fit to the 
data within the —3 < r) < 3 region using Eq. [T3] (see text for 
details). The shaded band represent 90% C.L. systematic er- 
rors. The open points were obtained by the tracklet analysis 
in the range \rj\ < 1. 



g(N part ) = l + 0129Nl / J t , (11) 

as shown by solid curves in Fig.[T4"h. Again, Fig. 114b dis- 
plays the ratio of data to this fit, which exhibits only a 
small deviation from unity. The line represents the aver- 
age value (R=1.0026)of this ratio over all data points and 
the shaded band the corresponding standard deviation of 
cr=0.047. 



B. Pseudorapidity distributions 

The final pseudorapidity distributions, dNch/di], are 
obtained by a simple equal weight averaging of the results 
from the hit-counting and energy deposition methods de- 
scribed in Sect. IIII Gl and Sect. IIIIH1 respectively. The 
distributions are shown as solid points in Figs. [TSIITSI for 



Au+Au collisions and in Figs. IT9"ll2"Tl for Cu+Cu collisions 
for different centrality bins. The shaded bands represent 
the range of systematic errors to the 90% confidence- 
level. 

All dN c h/dr] distributions exhibit a plateau around 
r\ ~ 0, the range of which increases with collision en- 
ergy followed by a smooth fall-off to higher values of \r]\. 
The fall-off, which is associated with the extended lon- 
gitudinal scaling region (see Sect. IIV El and Ref. [To}), is 
increasing in range with energy. It is also apparent that 
the level of the central plateau increases with both cen- 
trality and collision energy. Note that an earlier analysis 
of the 19.6 GeV Au+Au data gave up to 10% lower 
values. A re-analysis of this data set using an improved 
centrality determination and an improved dead-channel 
map led to the more reliable measurement presented here. 
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FIG. 16: dN c h/drj vs r\ (solid points) for eleven centrality 
bins representing 50% of the total cross section for Au + Au 
collisions at ^/s NN =62A GeV. The solid curve is a fit to the 
data within the —4.2 < r\ < 4.2 region using Eq. [12] (see text 
for details). The shaded band represent 90% C.L. systematic 
errors. The open points were obtained by the tracklet analysis 
in the range \r/\ < 1. 



1. Total charged particle multiplicity 

Nearly all the charged particles fall within the accep- 
tance of the PHOBOS multiplicity array. Therefore, it 
is appropriate to integrate the dNch/dr] distributions in 
order to estimate the total number of charged particles 
emitted in the collision, which was done using three dif- 
ferent methods as illustrated in FigfJH The first esti- 
mate, iVc/i||?7|<5.4) is a simple integration of dN cn /dri over 
the rj acceptance, — 5.4 < 77 < 5.4 of the multiplicity ar- 
ray corresponding to the shaded areas in Fig. [511 The re- 
sults of this analysis are listed in column 4 of Tables IVIII 
and IVIIII (see Appendix) and plotted as open circles in 
Figs. [Ml and [Ml 

By inspection of Figs. [15] and Q15] showing the 
dN c h/dri distributions for Au+Au and Cu+Cu collisions 
at y/T^=19.6 and 22.4 GeV, respectively, it is apparent, 



however, that large high-?7 tails develop as one moves to- 
ward more peripheral collisions. One possibility is that 
these tails represent charged particles emitted from colli- 
sion spectators, which travel in a very forward direction 
after being sheared off from the participant part of the 
incoming nuclei during the collision. The fact that the 
tails become more prominent for peripheral collisions and 
lower collision energies is consistent with this picture. 
Under this assumption, we estimate the multiplicity of 
charged particles originating from the collision zone by 
excluding these tails using the following procedure. We 
have found that for central collisions of Au+Au at 62.4, 
130 and 200 GeV and Cu+Cu at 62.4 and 200 GeV, the 
dN c h I drj distributions are well reproduced by 



dN ch _ cy/1 - l/(a cosh 77)2 
drj ~ l + e (H-' 3 )/ Q ' 1 ' 

where c, a, /?, and a are fit parameters. The quality of 
such fits, shown as solid curves, can be examined in e.g. 
Figs. Hlk and Mk- 

It has been demonstrated 0, EH that Au+Au and 
Cu+Cu collisions rather accurately obey extended lon- 
gitudinal scaling (a.k.a. limiting fragmentation scaling), 
even for non-central collisions, such that the multiplicity 
in the range around 77 = ybeam is independent of colli- 
sion energy. Consequently, we expect that the particle 
production in the range -ybeam < r) < ybeam is strongly 
dominated by participant collisions (i.e. excluding spec- 
tator emission). A fit to the data within this pseudo- 
rapidity range at all collision energies was therefore as- 
sumed to describe the contribution to the total charged 
particle multiplicity from participant collisions also out- 
side this interval. The solid black curves in Figs. Qj)] - 
[2T1 represent these fits to the data using the functional 
form and fit interval specified in the figure captions. The 
estimated total charged particle multiplicities for partic- 
ipant collisions obtained from this procedure are denoted 
N^ h and are listed in Tables IVlTl and [Villi (see Appendix) 
and shown as solid circles in Figs. |2"31 and |2"41 The quan- 
tity N^ h should represent the most reliable estimate of 
the total number of charged particles emitted from the 
overlap zone between the colliding nuclei. 

The third estimate, Nfj t represents an extrapolation 
outside the measured r\ region, which does not exclude 
contributions from spectator emission. This method was 
used in all earlier PHOBOS publications, e.g. @, fTTl — 
\\a\ . For the lowest collision energy, 19.6 GeV and 22.4 
GeV for Au + Au and Cu + Cu collisions, respectively, an 
average of iVc/i||?7|<5.4 an d N^ h is used. This corresponds 
to the area underneath the dashed curve in Fig. [22k . 
For higher collision energies, the experimental data (solid 
points in Fig. I2"2"b.d) are extended outside the measured 
77 region by shifting the low energy distributions by A77 = 
±[Ubea m - Vo], where ybeam is the beam rapidity at the 
collision energy in question and yo is the beam rapidity 
at the lowest collision energy. This extension of the data 
is thus based on the assumption that dN c h/dr] for rj > 
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FIG. 17: Same as Fig.[l6]but for Au+Au collisions at ^/s NN — 
130 GeV. The solid curves represent best fits to the data over 
the region -4.9< rj <4.9 using Eq.[l2](see text) and the shaded 
regions represent the systematic error band at 90% confidence 
limit. The open points were obtained by the tracklet analysis 
in the range \rj\ < 1. 



FIG. 18: Same as Fig. I17l but for Au + Au collisions at ^/s NJV 
= 200 GeV. The solid curves represent best fits to the data 
over the full 77 range using Eq. [12] and the shaded regions 
represent 90% C.L. systematic errors. The open points were 
obtained by the tracklet analysis in the range j?7| < 1. 



Vbeam is independent of collision energy - the limiting 
fragmentation hypothesis. The N^f estimate is thus the 
integral of the average of this extended distribution and 
the fit using Eq. [T27 dashed curve in Fig. I2"2"b). 

Tables IVlTl and [VI 1 1 1 ( see Appendix) summarize the to- 
tal charged particle multiplicity results for Au+Au and 
Cu+Cu collisions, respectively. The estimated average 
number of participants associated with each centrality 
bin was obtained from Glauber model (Monte Carlo ver- 
sion) [2^ and listed in column two. Column three lists 
the full width at half maximum (FWHM) of the dN ch /dri 
distributions, whereas the three total multiplicity esti- 
mates discussed above are listed in columns 4-6. 

The upper panels of Figs. [23] and [M]display the values 
N% h (solid points) and N c h\\ v \<5A (open circles). In all 
cases one observes participant scaling (Ref. [28j), an es- 
sentially linear dependence on (N part }. This is illustrated 
more clearly in the middle panels, where the participant- 



scaled results, dN c h/dr]/(N par t/2} are seen to be essen- 
tially independent of (N part /2} and exceeding the values 
obtained in pp/pp collisions. We observe that this quan- 
tity is almost constant with collision centrality. It is in- 
teresting to note that the normalized particle production 
in heavy-ion collisions is larger by about 40% than those 
of pp collisions (solid squares) [43| and pp collisions (solid 
diamonds) [42|. 

The widths of the dA^/d^-distributions, represented 
by the Full Width at Half Maximum (FWHM) are shown 
in the bottom panels in Figs.[23]and [M]for Au + Au and 
Cu + Cu collisions, respectively. The FWHM exhibit a 
decline with centrality, which indicates that the increased 
particle production with centrality preferentially occurs 
in the midrapidity region. Note also that the FWHM 
for pp reactions at 200 GeV and pp reactions at 62.4 
GeV and 19.6 GeV follow the trend of the Au + Au data 
extrapolated to (N part ) = 2. A similar trend is found for 
Cu + Cu collisions. 
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FIG. 19: dN c h/dr] vs 77 (solid points) for twelve centrality 
bins representing 55% of the total cross section for Cu + Cu 
collisions at % /s JVJV =22.4 GeV. The solid curve is a symmet- 
ric double quasi- Gaussian function fit to the data within the 
—3.2 < 77 < 3.2 region. The shaded band represent 90% 
C.L. systematic errors. The open points were obtained by 
the tracklet analysis in the range \rj\ < 1. 



FIG. 20: dN c h/dri vs 77 (solid points) for twelve centrality 
bins representing 55% of the total cross section for Cu + Cu 
collisions at ^js NN =62A GeV. The solid curve is a fit to the 
data within the —4.2 < 77 < 4.2 region using Eq. [T2] The 
shaded band represent 90% C.L. systematic errors. The open 
points were obtained by the tracklet analysis in the range 

H<i. 



2. Energy dependence 

In this sub-section we describe a simple semi-empirical 
expression for the total charged particle multiplicity in 
central heavy-ion collisions. It is motivated by the ob- 
servation that these distributions are largely character- 
ized by a midrapidity plateau, the height of which is 
well described by the empirical relation Eq. [7] and the 
fact that the width increases with collision energy such 
that the extended longitudinal scaling [7[ is fulfilled. A 
linear fit to the Au + Au data at 0-3% centrality in 
the region 2.0 < rj < 5.0 shows that it is well repro- 
duced by the relation dN c h/dr) = a(yb eam +f]o — rj), and 
dN ch /drj = ad/beam+m+v) f° r -5-0 < V < -2.0, where 
?/o=0.11 represents a small pseudorapidity offset and the 
slope in the fragmentation region is a=205. Since the 
mid-rapidity region is approximately flat, it is appropri- 
ate to use a trapezoidal shape of the dN/drj distributions 



which leads to an expression of the form 



NT 



dN, 



ch\Q 



dl] 



(2?7 + 2y be 



part, 



dN, 



h - 



2a 



dr) 



(13) 



where the term 2yb e am accounts for the increased width 
of the distribution as a function of collision energy, see 
Ref. for a similar analysis of the 0-6% centrality bin. 
Using the approximation ybeam — ^ ms iviv — M^oC 2 ), 
which is valid for ^/s NN ^ mo (mo being the nucleon 
mass) one obtains 



NT 



(N part /2} 



0.26(ln s NN ) 2 + 0.01 In s NN -0.28, (14) 



which is compared to experimental data for Au+Au, Pb+ 
Pb, and Cu + Cu collisions in Fig. [55] (dashed curves). 
This expression reproduces the PHOBOS data quite well, 
as expected, but underestimates the charged particle pro- 
duction at lower energies, presumably due to a break- 
down of the trapeziodal shape approximation at these 
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FIG. 21: Same as Fig.[20]but for ^/s^7 = 200 GeV. The solid 
curve is a fit to the data over the full r\ region using Eq. [12] 
The open points were obtained by the tracklet analysis for 
events in the range |?7| < 1. 



lower energies. This derivation does, however, predict 
the leading (lns N]v ) 2 term. We find that the expression 
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FIG. 22: Illustration of the methods used for estimating the 
total number of charged particles from the measured dN c h /dn 
vs T) (solid points) distributions. The solid curve represents a 
fit to the data within the region —yteam < rj < ybeam using 
the functional form of Eq. ll2l The open circles and the shaded 
regions are explained in the text. 



as discussed above, no centrality dependence outside of 
error bars is observed. The degree to which this quan- 
tity depends only on the collision energy is illustrated 
in Fig. [26l where N^ h / (N part /2) (solid points with er- 
ror ellipses) and N c h\\ n \<5A/{Npart/2) (open circles) are 
compared to the predictions obtained by Eq. [TS] (solid 
lines). The lower panels shows the ratio between the 
data, N^J {Np ar t/2), and the fit (Ea.fl5]l. Here the solid 
lines are the average of all points and the grey band repre- 
sents the standard deviation, which amounts to er= 0.036 
and g— 0.057 for Au + Au and Cu + Cu collisions, re- 
spectively. 



d+Au collisions 



NX 



(N par t/2} 



0.26(ln SjvJV ) 2 + 0.12, 



(15) 



gives an excellent description of the overall charged parti- 
cle production in heavy-ion collisions over the full energy 
range from y/s NN = 2.4 GeV to ^/s NN =200 GeV (solid 
lines in Fig. 1251) . 



3. Factorization 

Similar to what was shown for the midrapidity mul- 
tiplicity (Sect. IV. A. 2) the total participant-scaled 
charged particle multiplicity, dN c } l /drj/(N part /2), ap- 
pears to exhibit factorization of the centrality and en- 
ergy dependence, albeit of a somewhat trivial sort since, 



The primary reason for measuring d + Au collisions in 
RHIC Run 3 was to obtain high px spectra for a com- 
parison to Au + Au collisions and literature data on pp 
collisions to determine whether the observed high px sup- 
pression in Au + Au collisions is an entrance channel ef- 
fect (gluon saturation) or associated with high energy 
loss rates of partons traversing a color-charged medium. 
Concurrent with these studies, PHOBOS recorded the 
charged-particle multiplicity, the results of which are 
shown in Fig. 1271 for five centrality bins and a minimum- 
bias trigger. One observes that the dN c h/drj distribu- 
tion is strongly asymmetric for the most central collisions 
(panel e) but approaches symmetry for the most periph- 
eral collisions (panel a) . The asymmetry is expected from 
momentum conservation, since a larger number of partic- 
ipants are associated with the incoming Au nucleus trav- 
eling in the negative r\ direction, especially for central 
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FIG. 23: The total charged particle multiplicity, N c h, 
the total charged particle multiplicity per participant pair, 
2N c h/ {Npart} and the full width at half maximum, FWHM 
are shown as a function of centrality for each energy for 
Au + Au collisions. The first two quantities are shown for 
both the measured region A r c h|| r| | <5 . 4 (open points) and based 
on extrapolations, N% h (solid circles), see text for details. Cor- 
responding data for inelastic pp collisions [43j] at 200 GeV 
and pp collisions !42| are shown as solid squares, whereas the 
crosses are from the present work and solid diamonds rep- 
resent an overall fit N ch = —0.42 + 4.69s 0155 for inelastic 
collisions H3|. The FWHM for 200 GeV pp and 62.4 GeV 
pp collisions (solid squares) were obtained from Refs. [43|] and 
Ref. p2|. respectively, whereas the value (solid diamond) for 
19.6 GeV represents an extrapolation based on data in the 
latter work. No FWHM data are available for 130 GeV pp 
collisions. Error bars and ellipses represent systematic 90% 
C.L. errors. 
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FIG. 24: Same as Fig.[23]but for Cu +Cu collisions. The data 
for 22.4 GeV pp collisions (solid diamonds) were obtained as 
for the 19.6 GeV data shown in Fig. [23] 
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collisions, as predicted in Monte Carlo Glauber model 
calculations, see Table IIIII The total charged-particle 
multiplicity, N c h\\n\<5.ii within the -5.4< r\ <5.4 accep- 
tance region is also listed along with an estimate of the 
total charged particle multiplicity, iV*?*, which includes 
an estimate of the contribution from the unmeasured re- 
gion. As detailed in Ref. [l(| this is obtained from lower 
energy p+A collisions by assuming extended longitudinal 
scaling. Also discussed in Ref. [10| (see Fig. 4a) is the 
observation that the participant-scaled total multiplicity, 
Ntff / (Np art /2), for d + Au collisions is commensurate 
with that observed for pp collisions, see e.g. Fig. 13 of 
Ref. The additional (~40%) enhancement seen for 

Au+Au and Cu+Cu systems does not appear for d+Au. 



FIG. 25: Doubly logarithmic (panel a) and linear (panel b) 
plots of N v ch /(N part /2) vs. (Ins) 2 . The PHOBOS results are 
for 0-3% central collisions. The trend is extended to lower 
energies with data from the SPS 22;j and AGS The 
dashed and solid curves were calculated using Eqs, I14l and ll5l 
respectively. 



D. p+p collisions 

Inelastic proton-proton collisions were measured at two 
energies, 200 GeV and 410 GeV. The 200 GeV data were 
measured, in part, in order to provide a baseline against 
which to identify and study the special effects associated 
with heavy-ion collisions in observables such as charged 
particle multiplicities, collective flow, and high px sup- 
pression in particle spectra. In addition, pp data were 
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FIG. 26: Two estimates of the total charged particle 
multiplicity per participant pair, 2N ch \^ <5A / {N pa rt} (open 
points) and 2N^ h / {N par t} (solid points), are shown as a func- 
tion of centrality for each energy for Au + Au (panel a) and 
Cu + Cu (panel c) collisions. Horizontal lines represent cen- 
trality independent levels predicted by Eq. II 5 1 The lower 
panels show the ratio of data to the Eq. [TS] prediction. The 
solid line is the average of all data points and the grey-shaded 
area the la standard deviation. Corresponding data for pp/pp 
collisions are shown as stars, see Figs. [23] and [24] for details. 



TABLE III: Summary of the Monte Carlo Glauber model 
predictions of the number of participating nucleons, (N^Xrt) 
and (Np art ) associated with the incoming Au and deuteron 
nuclei, respectively. The total number of charged particles 



emitted within the PHOBOS acceptance, N c h 



M<5-4 



and an 



estimate of the total multiplicity including the unmeasured 
region, iV*£ 4 (see text for details) are also listed. Errors are 
90% C.L. systematic errors. 



Cent. (%) 


(AjAu \ 
\^ v part / 


impart) N c h\\ v \<5.4, 




0-20 


13.5±1.0 


2.0±0.1 


157±10 


167+it 


20-40 


8.9±0.7 


1.9±0.1 


109±7 


115+8° 


40-60 


5.4±0.6 


1.7±0.2 


74±5 


7711 


60-80 


2.9±0.5 


1.4±0.2 


46±3 


48+j 


80-100 


1.6±0.4 


1.1±0.2 


28±3 


29 + 3 


Min-Bias 


6.6±0.5 


1.7±0.1 


82±6 


87t^ 



collected at the highest RHIC energy, namely s/s = 410 
GeV. 

Only the hit-counting method was applied in the anal- 
ysis of data for pp collisions. Using HIJING Monte Carlo 
simulations, the particle yields were corrected for detec- 
tor acceptance, secondary particles produced in the ma- 
terial of the detector, and particles resulting from weak 
decays. These corrections have similar dependences on r\ 
as those found for heavy ion collisions. Pscudorapidity 
distributions were extracted in bins of the total observed 
charged-particle multiplicity. The distribution found for 
each bin was corrected for the effects of triggering and 
vertexing efficiency, which were both found to be strongly 
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FIG. 27: The charged particle multiplicity, dN c h/dri, is shown 
for five centrality bins for d + Au collisions (panels a-e). The 
grey-shaded area represents the la standard deviation. The 
minimum-bias distribution is shown in panel f. The pseudora- 
pidity, rj, is given relative to the deuteron beam, which travels 
in the positive z direction. 



multiplicity dependent. The analysis compared the num- 
ber of input Monte Carlo events to the number surviving 
all event criteria so separate triggering and vertex effi- 
ciencies were not extracted. Individually corrected distri- 
butions were then combined, weighted by the efficiency 
corrected number of events in each multiplicity bin, to 
generate the final average dN/dry. Systematic uncertain- 
ties, which were much larger than the statistical ones, 
were found by varying the cuts used to select events and 
hits, by comparing the results for positive and negative 
pseudorapidity, and by considering the differences found 
between the hit-counting and other techniques when ap- 
plied to heavy ion data. Systematic uncertainties due to 
all sources were then added in quadrature. 

The charged-particle multiplicity distributions in pseu- 
dorapidity are presented in Fig. [25] for inelastic pp colli- 
sions [26| . The grey-shaded area around the data points 
shows the 90% confidence limit systematic error. The 
midrapidity and total charged particle multiplicity for 
minimum bias inelastic pp collisions, extracted using the 
symmetric double-Gaussian fits to the data shown as 
solid curves in Fig. [25] are listed in Table IIV1 



E. Extended longitudinal scaling 

It is a well-known phenomenon that, at sufficiently 
high energy, the particle production in the rapidity re- 
gion of either collision partner becomes largely indepen- 
dent of the collision energy. This effect is referred to as 
"limiting fragmentation scaling" [52j and it has been ob- 
served in p + p,p + A and heavy-ion collisions [7]. The 
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TABLE IV: Table of the midrapidity charged-particle den- 
sity dN c h/dr], the total number of charged particles emitted 
within the PHOBOS acceptance, N c h\\m\<sA; m minimum- 
bias, inelastic p + p collisions. Also listed is an estimate of the 
total multiplicity including the unmeasured region, iV*jf us- 
ing a three-parameter double-Gaussian fit to the data. Errors 
are 90% C.L. systematic errors. 
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FIG. 28: The charged particle multiplicity, dN c h/dr], is shown 
for 200 GeV (panel a) and 410 GeV (panel b) for pp inelastic 
collisions. The grey-shaded area represents the 90% confi- 
dence limit systematic error. The solid curves are symmetric 
double- Gaussian fits to the data used to derive the total mul- 
tiplcities listed in Table IIVI 



term "extended longitudinal scaling" has also been used 
to describe this phenomenon, since this scaling feature 
covers a more extended region of 77 than expected from 
the hypothesis of limiting fragmentation and because it 
also appears to apply to other observables in heavy-ion 
collisions, e.g. the magnitude of elliptic flow (53[. Al- 
though the original concept refers to the particle produc- 
tion as a function of rapidity, y, it can be shown that the 
scaling also holds for dNch/dr]. This results from the fact 
that for 1 77 1 3> 1 the relation pt sinh 77 = mrsinhy leads 
to y ~ ln(px /tut) + rj, where the term ln(pT /fir) leads 
to a very small correction that is not taken into account. 
Similarly, the Jacobian px cosh 77/ y m 2 + px cosh 2 77 as- 
sociated with the transformation dN c h j dy to dN c h /dr] is 
close to unity in the 77-region of interest and may safely 
be ignored. In Figure UHl the extent to which the lim- 
iting fragmentation scaling is valid for Au + Au (panels 
a,b), Cu + Cu (panels c,d) and p + p collisions (panels 
e,f) is shown in both the positive and negative 77 regions. 
Further discussion of "limiting fragmentation scaling" in 




-2 2 4 
^ — ^ ^beam 



6-6 -4 -2 
^ — ^ ^beam 



FIG. 29: Illustration of the extended longitudinal scaling 
(a.k.a. limiting fragmentation scaling) for Au + Au (pan- 
els a, b), Cu + Cu (panels c, d), and p + p collisions (panels 
e, f). The systematic errors (90% C.L.) are shown as shaded 
areas for the highest energy only for each system. The arrows 
indicate the location of mid- rapidity (77=0). 



heavy-ion collisions can be found in Refs. 0, [Hj]. 



V. CONCLUSIONS 

We have measured the multiplicity of charged parti- 
cles in Au + Au collisions at nucleon-nucleon center-of- 
mass energies of ys~^7=19.6, 56, 62.4, 130, and 200 GeV, 
Cu + Cu collisions at ^fs^=22A, 62.4, and 200 GeV, 
^/s JVJV =200 GeV d + Au collisions, and p + p collisions 
at A /s Niv =200 and 410 GeV at the Relativistic Heavy- 
Ion Collider at Brookhaven National Laboratory. For 
the Au + Au system, the complete distributions in pseu- 
dorapidity were measured at y/57^7^19.6, 62.4, 130, and 
200 GeV for 10-11 centrality bins corresponding to 45- 
50% most central collisions, whereas the measurement 
at 56 GeV yielded only the charged particle multiplicity 
at midrapidity for the 6% most central collisions. The 
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midrapidity measurements, based on the tracklet analy- 
sis method were extended up to 70% of the cross section 
in up to 15 centrality bins. For Cu + Cu, the pseudora- 
pidity distributions were measured at ^/s NN =22.4, 62.4, 
and 200 GeV in 12 centrality bins corresponding to 55% 
most central collisions. 

Measurements were also carried out for the asymmetric 
d + Au system at ^§"^"=200 GeV, and the multiplicity 
distributions are reported for five centrality bins as well 
as minimum bias events. Also elementary p + p collisions 
were measured at t/s=200 and 410 GeV and the results 
given in the present work. 

The multiplicity results were in all cases derived from 
three different methods of analysis, which all agree within 
systematic errors. The multiplicities were all corrected 
for missing low px particles and the contribution from 
weak decays have been subtracted. At midrapidity, we 
observe a smooth decrease of charged particle production 
per participant nucleon pair when going to more periph- 
eral collisions trending towards the results from pp and pp 
collisions for the most peripheral collisions at the same 
energies. Surprisingly, the midrapidity production can 
be factorized into separate energy and centrality depen- 
dences for both the Au + Au and Cu + Cu data. Com- 
bined with lower energy data from the SPS (CERN) and 
AGS (Brookhaven) we find an approximate logarithmic 
energy dependence of the midrapidity charged particle 
production. 

A unique feature of the PHOBOS experiment is that 
the pseudorapidity distributions were measured over a 
very wide range of —5.4 < r\ < 5.4. Three estimates of 
the total charged particle multiplicity in Au + Au and 
Cu + Cu collisions were performed. One method at- 
tempts to exclude the contribution from excited spec- 
tators, which appear to contribute substantially at the 
highest pseudorapidities for the lowest collision energies, 
by using an analytical function fit to the data within the 
pseudorapidity range -y beam < rj < ybeam- By extrapo- 
lation into the unmeasured 77-region, of the order 1-5% in 
total multiplicity, it is thus possible to estimate the total 
charged particle multiplicity generated in participant col- 
lisions. When scaled to the number of participant pairs, 
this multiplicity is constant as a function of centrality 
and about 40% larger than those observed for pp and 
pp collisions. Approximating the central pseudorapidity 
distributions by a trapezoidal shape, the height of which 
is given by the logarithmic midrapidity dependence, and 
the width is derived from the fact that the wings of the 
distribution follow the extended longitudinal scaling, we 
have found that the total charged particle production 
follows a simple (Ins) 2 scaling, which applies over the 
full energy range from the lowest AGS energies y/s=2A 
GeV to the most energetic RHIC collisions v / s=200 GeV. 
However, an interesting and statistically significant de- 
parture from the otherwise smooth energy dependence 
of the total multiplicity and midrapidity density for all 
centralities is seen for ^/s NJV ~60 GeV, which may be an 
indication of curvature. 



The widths of the pseudorapidity distributions are 
found to decrease with centrality at a rate that is approx- 
imately compensated for by the increase in the midrapid- 
ity plateau with centrality such that the total multiplicity 
is found to be almost constant when normalized to the 
number of participant pairs (N part )/2. 

Finally, it should be noted that most of the prominent 
trends seen in the PHOBOS multiplicity data are not spe- 
cific to A + A collisions at RHIC energies. For example, 
Np ar t scaling is seen in all hadron + A collisions [ID, [28[ , 
extended longitudinal scaling is seen in p + p and e + e~ 
collisions [l4| , and also in p + A collisions [l8|, H3| , and 
the logarithmic rise of the midrapidity particle density is 
seen in A + A collisions at lower energies (see Fig. [T2")) . 

A complete tabulation of the pseudorapidity distribu- 
tion data shown in Figs. [T51I2T1 |2"T1 and [35] as well as 
various derived quantities not listed in this article, may 
be obtained from the Electronic Physics Auxiliary Publi- 
cation Service of the American Institute of Physics. See 
EPAPS Document No. [number will be inserted by pub- 
lisher ]. A direct link to this document may be found 
in the online article's HTML reference section. The 
document may also be reached via the EPAPS home- 
page (http:/ /www. aip.org/pubservs/epaps. html) or from 
ftp.aip.org in the directory /epaps/. See the EPAPS 
homepage for more information. 
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Appendix: Au+Au and Cu+Cu multiplicity tables 

In this appendix we give the tables of relevant param- 
eters for the PHOBOS charged particle multiplicity data 
for Au + Au and Cu + Cu collisions. The corresponding 
parameters for d + Au and p + p collisions are given in 
tables ITTil and ITVl in the main text. 
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TABLE V: Summary of the midrapidity Im<i charged 
particle multiplicity for Au+Au collisions obtained from the 
tracklet analysis. The data are listed as a function of central- 
ity expressed in percentage of the total reaction cross section 
for all four energies. Columns 2 and 4 list derived quantities, 
namely number of participants as well as the midrapidity den- 
sity normalized to the number of participant pairs N part /2. 
The errors are systematic errors at 90% C. L.; statistical er- 
rors are negligible. Note, the table continues overleaf. 



Au+Au 


N N 


19.6 GeV 


Vbeam— 3.036 


Bin 


N part 


dJV | 

dr/ 1 1 V 1 < 1 


dN ch/ d V\ |,,<1 
1 Y part I ^ 


0-3% 


351+11 


331+24 


1 89+0 15 


3-6% 


322+10 


297+22 


1.84+0.15 


6-10% 


286+9 


260+20 


1.82+0.15 


10-15% 


247+8 


216+16 


1.76+0.14 


15-20% 


206+8 


181+14 


1 75+0 15 


20-25% 


171+7 


148+11 


1.73+0.15 


25-30% 


142+7 


121+9 


1.70+0.15 


30-35% 


117+7 


97+7 


1.65+0.16 


35-40% 


95+7 


78+6 


1.64+0.17 


40-45% 


74+6 


59+4 


1.61+0.18 


45-50% 


58+3 


45+3 


1.56+0.19 


50-55% 


45+3 


35+3 


1.55+0.20 


55-60% 


34+3 


26+2 


1.55+0.22 


Au+Au 






iiu —A 1 Qfi 


Bin 


7V r 

J V part 


dN 1 

dr, 1 1 ^7 1 < ! 


d ^ah/ dr l\\ri<l\ 
Npart / 2 


0-3% 


356+11 


492+36 


2.76+0.23 


3-6% 


325+10 


433+32 


2.67+0.22 


6-10% 


288+9 


377+28 


2.62+0.21 


10-15% 


248+8 


316+23 


2.55+0.21 


15-20% 


209+7 


260+19 


2.50+0.20 


20-25% 


174+7 


212+15 


2.44+0.21 


25-30% 


145+7 


174+13 


2.41+0.21 


30-35% 


119+7 


140+10 


2.35+0.22 


35-40% 


98+7 


111+8 


2.28+0.23 


40-45% 


78+6 


87+6 


2.24+0.25 


45-50% 


62+6 


67+5 


2.16+0.26 


50-55% 


48+5 


50+4 


2.10+0.27 


50-60% 


36+4 


36+3 


2.01+0.28 


60-65% 


27+3 


25+2 


1.91+0.28 


65-70% 


19+3 


17+1 


1.77+0.27 
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TABLE V: Continued. TABLE VI: Same as Table El but for Cu+Cu collisions. 



Au+Au y/s^= 


=130 GeV 


y(, C am=4.930 


Cu+Cu 




=22.4 GeV 


yf>eam=3.170 


.Dill 


7V i 

1 'part 


Ml, , , 


dN ch/ d v\ \ v <i\ 


Bin 


AT 


dJV 


dJV ch/ d ll|T7<l 


d-q M<1 




d?j IM<1 


N va rt 1 2 


0-3% 


355±12 


613±24 


O A f 1 r\ 1 1 

3.45+0.17 


0-3% 


103+ 3 


101+8 


1.96+0.16 


3-6% 
u- ±u /o 


330±10 

one: i n 


545±21 
472+18 


3.31±0.16 
3.20+0.16 


3-6% 
6-10% 


95+ 3 
86+ 3 


90+7 
80+6 


1.90+0.15 
1.86+0.15 


10-15% 


254±8 


QQQ-|--| c: 




10-1570 


74+ 3 


68+5 


l.83±0.lo 


15-20% 


214±8 


327±13 


3.06±0.16 


10- zU /o 


Doit o 


t ±4 


l.oZztU. 10 


20-25% 


179±7 


m A 1 -1 I 

Z 1 4zhll 


Q Ofi_l_n 1 7 


20-25% 


53+ 3 


48+4 


1.81+0.17 


25-30% 


148±6 


220±8 


2.96±0.17 


25-30% 


44+ 3 


40+3 


1.80+0.19 


30-35% 


122±6 


180±7 


2.94±0.18 


30-35% 


37+ 3 


33+2 


1.78+0.21 


35-40% 
40-45% 


100±5 
80±5 


140±5 
110±4 


2.80±0.18 
2.75±0.20 


35-40% 
40-45% 


30+ 3 
24+ 3 


27+2 
22+2 


1.79+0.24 
1.76+0.26 


45-50% 


63±4 


83±3 


2.64±0.21 


45-50% 


20+ 3 


18+1 


1.77+0.29 


Au+Au % /s JVJV = 


=200 GeV 


Ubeam — 5. 361 


Cu+Cu 




=62.4 GeV 


yi,eam=4.196 


Bin 


A^ part 


dJV 1 


dA, ch./ d ')l|,,<l 


.Dill 


M i 

+ ' part 


dN \, 


dN ch/ d n\\n<l\ 


dq N»?I<1 


Npart / 2 


drj I|»7|<1 


AW/ 2 


0-3% 


361±11 


691±52 


3.82±0.31 


0-3% 


106+3 


138+10 


2.64+0.21 


3-6% 


331±10 


619±46 


3.74±0.30 


3-6% 


97+3 


123+ 9 


2.55+0.21 


6-10% 


297±9 


540±41 


3.64±0.30 


6-10% 


88+3 


108+ 8 


2.46+0.20 


10-15% 


255±8 


465±35 


3.65±0.30 


lO-lo/o 


/6±3 


no_L_ 1 
9z± 1 


o /in i a on 


15-20% 


215±7 


o o a I on 

384±29 


3.57+0.29 


lO-zU/o 


o rr I o 
00±O 


iia~ a 


o QTin on 
1.6 ( dzU.zU 


20-25% 


180±7 


313±24 


3.47±0.30 


20-25% 


55+3 


64+ 5 


2.33+0.22 


25-30% 


150±6 


257±19 


3.42±0.29 


25-30% 


47+3 


52+ 4 


2.25+0.23 


30-35% 


124±6 


20o±16 


o 97±n on 


30-35% 


38+3 


43+ 3 


2.26+0.26 


35-40% 


101±6 


165±12 


3.25±0.31 


35-40% 


32+3 


35+ 3 


2.22+0.28 


40-45% 


82±6 


133±10 


3.25±0.34 


40-45% 


26+3 


28+ 2 


2.21+0.32 


45-50% 


65±6 


100±8 


3.10±0.38 


45-50% 


21+3 


23+ 2 


2.21+0.35 


50-55% 


49±5 


73±5 


2.98±0.37 


Cu+Cu 




=200 GeV 


ybeam=5.361 


55-60% 


37±4 


54±4 


2.88±0.39 


.Dill 


M 


dJV 1, 

dr, IM<1 


dN ch/ d n\\n<l\ 


60-65% 
65-70% 


28±3 
20±3 


38±3 
27±2 


2.78±0.40 
2.68±0.41 


impart 


N va rt 1 2 


0-3% 
3-6% 


108+4 
101+3 


198+15 
175+13 


3.66+0.29 
3.48+0.28 


















6-10% 


91+3 


155+12 


3.42+0.28 










i n i ccoy 
lU-lo/o 


id I o 

/9±o 




9 QQ_l_n 07 

o.ooztU.Z / 










15-20% 


67+3 


109+ 8 


3.26+0.28 










20-25% 


57+3 


91+ 7 


3.21+0.29 










25-30% 


48+3 


75+ 6 


3.17+0.32 










30-35% 


40+3 


62+ 5 


3.15+0.35 










35-40% 


33+3 


50+ 4 


3.07+0.38 










40-45% 


27+3 


40+ 3 


3.04+0.43 










45-50% 


22+3 


32+ 2 


2.97+0.46 
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TABLE VII: Summary of the total charged particle multiplic- 
ity estimates for Au+Au collisions obtained in three different 
ways. iVcfc||„i<5.4 denotes the multiplicity observed within 
the acceptance region \r]\ < 5.4 obtained from the single layer 
analysis, whereas NT and N* h represent extrapolations into 
the unmeasured r\ region, either including or excluding con- 
tributions from spectator emission. See text for details. The 
data are listed as a function of centrality expressed in per- 
centage of the total reaction cross section for all four energies. 
Also listed are the full width at half maximum, FWHM, of 
the dN/drj distributions and the derived quantity N part . Only 
systematic errors at 90% C. L. are given since statistical errors 
are negligible. 



TABLE VIII: Same as Table |VlI] but for Cu+Cu collisions 



7=19.6 GeV 



Vbeam 



=3.036 



Bin 


■^V part 


FWHM 


N c h\ |t;|<5.4 


Kk 


Kh 


0-3% 


353±11 


4.80t^ 


1682±115 


1676±115 


1669±115 


3-6% 


323±10 


4.86ig;» 


1531±111 


1522±111 


1512±111 


6-10% 


286±9 


4.94t°j 


1367±97 


1352±97 


1337±97 


10-15% 


246±8 


5.03+oi 


1182±86 


1162±86 


1145±86 


15-20% 


206±8 


5.101$ 


1014±78 


989±78 


973±78 


20-25% 


172±7 


5.19ir 3 


866±71 


836±71 


813±71 


25-30% 


142±7 


O.Zo Q 3 


735±68 


707±68 


684±68 


30-35% 


117±7 


5.37l° ;I 


617±66 


586±66 


563±66 


35-40% 


95±7 


5.411^ 


516±63 


482±63 


453±63 


40-45% 




418±60 


387±60 


362±60 





Au+Au yfs~^- 


=62.4 GeV 


2/teom=4.196 


Bin 


N part FWHM 


Nch | 1 17 1 < 5 . 4 


NT Kh 



0-3% 


349±11 


6A7±»i 


2935+147 


2988+149 


2971+149 


3-6% 


323±10 


6.49l°i 


2733+137 


2775+138 


2762+138 


6-10% 


288±8 


6.54l°i 


2448+122 


2489+124 


2471+124 


10-15% 


248±7 


6.62±g| 


2077+103 


2120+106 


2094+106 


15-20% 


209±7 


6.68lg"i 


1739+87 


1777+88 


1747+88 


20-25% 


174±7 




1448+72 


1485+74 


1451+74 


25-30% 


145±7 


6.831H 


1200+60 


1236+61 


1203+61 


30-35% 


120±7 


6.91+°! 


984+49 


1027+51 


986+51 


35-40% 


98±6 


6.961^ 


797+40 


840+42 


796+42 


40-45% 


79±6 


7.05t°;« 


644+32 


679+33 


641+33 


45-50% 


62±6 


7M±ll 


504+25 


532+26 


505+26 




Au+Au ^/s NN : 


=130 GeV 


Vbeam =4.930 


Bin 


N part 


FWHM 


N ch 1 1 77 1 < 5-4 


NT 


Kh 



0-3% 


354+12 


' - i0 -0.3 


4286+214 


4376+219 


4346+219 


3-6% 


327+10 


7 O1+0.3 
' ' zl -0.3 


3915+196 


4015+201 


3971+201 


6-10% 


298+9 


7.29±°;1 


3546+182 


3649+182 


3598+182 


10-15% 


258+8 


ry qq + 0.3 

' - oo -0.3 


2982+149 


3090+155 


3032+155 


15-20% 


215+8 


7-39i°;| 


2482+124 


2586+129 


2525+129 


20-25% 


178+7 




2056+103 


2164+108 


2096+108 


25-30% 


148+6 


7.52i°;l 


1686+84 


1793+90 


1719+90 


30-35% 


124+6 


7.60l°1 


1395+70 


1502+75 


1422+75 


35-40% 


100+5 


7.65ir 4 


1116+56 


1222+61 


1142+61 


40-45% 


80+5 


7.74±g;« 


885+44 


975+49 


908+49 


45-50% 


64+4 


7.81±H 


705+35 


782+39 


726+39 




Au+Au 


=200 GeV 


2/i>eam=5.361 


Bin 


N part 


FWHM 


Nch |r;|<5.4 


NT 





0-3% 
3-6% 
6-10% 
10-15% 
15-20% 
20-25% 
25-30% 
30-35% 
35-40% 
40-45% 
45-50% 



0.4 

+0.3 
0.4 

7 7 ,+0.4 

7.77l°; 4 4 



359+11 7.63 
330+10 7.64 
297+9 
256+8 
215+7 
181+7 
149+6 
123+6 
101+6 

82+6 

65+6 



7.83Z 
7.90^ 



0.4 
0.4 
+0.4 
0.4 

7.99±°1 
7.961H 
8.041^ 

8-ii±8:l 

8.231 ,' 5 , 



5159+258 
4753+238 
4198+210 
3598+180 
2999+150 
2496+125 
2057+103 
1685+84 
1354+68 
1096+55 
847+42 



5290+264 
4895+245 
4341+217 
3763+188 
3153+158 
2645+132 
2184+109 
1819+91 
1486+74 
1204+60 
951+48 



5261+264 
4854+245 
4288+217 
3687+188 
3075+158 
2564+132 
2119+109 
1742+91 
1403+74 
1137+60 
880+48 





Cu+Cu y/S^= 


22.4 GeV 






Bin 


part 


FWHM 


1 V fh 1 U 1 K A 

cn, | | rj \ <. o . 4 


NT 

ch 


ch 


0-3% 


103+3 


5 96 + i'* 


538+22 


535+23 


534+23 


3-6% 


95+3 


5 82 + !M* 

— 0.40 
" J,uo — 0.36 


485+20 


482+21 


480+21 


6-10% 


86+3 


435+18 


431+19 


429+19 


10-15% 


74+3 


. OO n <"io 
— V.ZQ 


380+16 


375+18 


372+18 


15-20% 


63+3 


r 4 o+0.40 
°-*°-0.30 


326+14 


320+15 


316+15 


20-25% 


53+3 


tr, 4 o+0.38 


279+12 


273+14 


268+14 


25-30% 


44+3 




237+10 


230+13 


226+13 


30-35% 


37+3 


c 9 o+0.40 
J - z °-0.26 


201+ 9 


194+12 


191+12 


35-40% 


30+3 


°' zu -0.24 


169+ 8 


162+12 


160+12 


40-45% 


25+3 


r i 9+0.34 
°' ±z -0.22 


141+ 7 


135+11 


133+11 


45-50% 


20+3 


5 n6 +0 - 30 

O.UD_g 24 


118+ 6 


112+11 


110+11 


50-55% 


16+3 


4 Q4+0-28 
^• yq: -0.22 


96+ 6 


92+11 


91+11 




CU + CU yfS~^= 


62.4 GeV 


Vb, 


zam 4. 196 


Bin 


part 


FWHM 


N c h\\ v \<5.4 


NT 


Kh 



0-3% 

3-6% 

6-10% 

10-15% 

15-20% 

20-25% 

25-30% 

30-35% 

35-40% 

40-45% 

45-50% 

50-55% 



102+3 6.46 
95+3 
88+3 
76+3 
65+3 
55+3 
47+3 
38+3 
32+3 
25+3 
21+3 
16+3 



0.55 

6 44+ 52 

fi 4R+0.54 
°- 4D -0.44 

6.501 048 
6.54" 1 



0.44 
+0.56 
0.48 

6.58+ 057 



0.50 
+0.58 
0.50 
+0.60 
-0.50 
+0.64 
0.50 
+0.84 
0.66 

6.80j:°;£ 

+0.86 
0.60 



6.62Z 
6.1 

6.701! 
6.76;! 



6.80 4 



824+36 833+36 834+36 

771+34 781+34 780+34 

710+31 721+32 717+32 

624+27 635+27 629+27 

530+24 541+24 534+24 

449+20 460+21 451+21 

375+17 386+17 377+17 

313+14 323+15 314+15 

261+12 270+13 261+13 

214+11 223+11 216+11 

174+ 9 183+ 9 175+ 9 

140+ 7 147+ 8 140+ 8 





Cu+C 


U \/ S NN- 


200 GeV 


J/i,eam=5.361 


Bin 


■A^ part 


FWHM 


^<a||jj|<5.4 


NT 


Kh 


0-3% 


106+3 


<-»U_0.48 


1506+67 


1541+70 


1542+70 


3-6% 


100+3 


7 Sfi+0- 54 


1370+66 


1407+68 


1407+68 


6-10% 


91+3 


7 QB + 0- 54 
'•°°-0.56 


1226+57 


1262+59 


1260+59 


10-15% 


79+3 


7.96l°:g 


1048+49 


1084+51 


1081+51 


15-20% 


67+3 


8.02±g;S 


882+42 


917+43 


912+43 


20-25% 


57+3 


8.04lg : g 


739+35 


771+38 


766+38 


25-30% 


47+3 


8.08l° : » 


616+30 


645+32 


640+32 


30-35% 


40+3 


8.i6±g : ?5 


512+25 


538+27 


533+27 


35-40% 


33+3 


s 99+0.70 
°- zz -0.72 


420+21 


444+23 


438+23 


40-45% 


27+3 


s 09+0. 70 
°- oz -0.82 


342+17 


364+19 


358+19 


45-50% 


22+3 


09+O.72 
°- oz -0.76 


274+14 


293+15 


288+15 


50-55% 


17+3 


a QO+0.80 

o-^-o.so 


218+12 


234+13 


230+13 



25 



[1] B. B. Back et al. (PHOBOS), Nucl. Instr. Meth. A499, 
603 (2003) 

[2] B. B. Back et al. (PHOBOS), Phys. Rev. Lett. 85, 3100 
(2000) 

[3] B. B. Back et al. (PHOBOS), Phys. Rev. Lett. 87, 102303 
(2001) 

[4] B. B. Back et al. (PHOBOS), Phys. Rev. C 65, 031901R 
(2002) 

[5] B. B. Back et al. (PHOBOS), Phys. Rev. Lett. 88, 022302 
(2002) 

[6] B. B. Back et al. (PHOBOS), Phys. Rev. C 65, 061901R 
(2001) 

[7] B. B. Back et al. (PHOBOS), Phys. Rev. Lett. 91, 052303 
(2003) 

[8] B. B. Back et al. (PHOBOS), Phys. Rev. Lett. 93, 082301 
(2004) 

[9] B. B. Back et al. (PHOBOS), Phys. Rev. C 70, 
021902(R) (2004) 
[10] B. B. Back et al. (PHOBOS), Phys. Rev. C 72,031901 
(2005) 

[11] B. B. Back et al. (PHOBOS), Phys. Rev. C 74,021901 
(2006) 

[12] B. B. Back et al. (PHOBOS), Phys. Rev. C 74,021902 
(2006) 

[13] B. Alver et al. (PHOBOS), Phys. Rev. Lett. 102,142301 
(2009) 

[14] B. B. Back et al. (PHOBOS), Nucl. Phys. A 757, 28 
(2005) 

[15] I. Arsene et al. (BRAHMS), Nucl.Phys. A757, 1 (2005) 
[16] Adler O, et al. (STAR), Phys. Rev. Lett., 87 112303 
(2001) 

[17] Adcox K., et al. (PHENIX), Phys. Rev. Lett. 86, 3500 
(2001) 

[18] J. E. Elias et al, [E178], Phys. Rev D22, 13 (1980) 
[19] D. H. Brick et al., [E565/570], Phys. Rev. D39, 2484 
(1989) 

[20] C DeMarzo et al., [NA5], Phys. Rev. D26, 1019 (1982) 

[21] S. Fredriksson et al., Phys. Rep. 144, 187 (1987) 

[22] S. V. Afanasiev et al., Phys. Rev. C 66, 054902 (2002); 

T. Anticic et al, Phys. Rev. C 69, 024902 (2004) 
[23] J.L.Klay et al. (E895), Phys. Rev. C68, 054905 (2003) 
[24] R. Bindel et al., Nucl. Instr. Meth. A 474, 38 (2001) 
[25] R.S. Hollis et al., Proceedings of Intl. Wkshp. Part. Mult. 

Rel. HI Coll. 1719 June 2004, Bari, Italy, J.Phys Conf. 

Ser. 5, 46 (2005) 
[26] J. F. Sagerer, Ph.D. Thesis, University of Illinois at 

Chicago (2008) 

[27] A. Bialas, M. Bleszynskiand, and W.Czyz, Nucl. Phys. 

Bill, 461 (1976) 
[28] W. Busza et al. (E178) Phys. Rev. Lett. 34, 836 (1975) 



[29] We use the HIJING Glauber model, which is a Monte 

Carlo method, to estimate the cross section integrals. 

This method is generally accepted to be more accurate 

than analytical methods. 
[30] R. Nouicer et al. (PHOBOS), Nucl. Instr. Meth. A461 

(2001)143. 

[31] E. Nygard et al, Nucl. Instr. Meth. A301(1991)506. 
[32] O. Toker et al., Nucl. Instr. Meth. A340(1994)572. 
[33] E. Griesmayer et al, IEEE Trans. Nucl. Sci. NS-48, 1565 
(2001) 

[34] K. Wozniak et al. (PHOBOS), Nucl. Instr. Meth. A566, 
185 (2006) 

[35] E. Garcia et al., Nucl. Instr. Meth. A570, 536 (2007) 
[36] X.-N. Wang and M. Gyulassy, Phys. Rev. D 44,3501 

(1992); 45, 844 (1992); M. Gyulassy and X.-N. Wang, 

Comp. Phys. Comm. 83, 307 (1994); X.-N. Wang, Phys. 

Rep. 280, 287 (1997); Nucl. Phys. A661, 609c (1999) 
[37] GEANT Detector Description and Simulation Tool v. 

3.21, CERN Program Library, W5013 (1994) 
[38] K. Werner, Phys. Rep. 232, 87 (1993) 
[39] H. Sorge, Phys. Rev. C 52, 3291 (1995) 
[40] M. Reuter, Ph. D. Thesis, University of Illinois at 

Chicago, (2001) 
[41] K. Hagiwara et al, Phys. Rev. D 66, 010001 (2002) 
[42] W. Thome et al. (ISR), Nucl. Phys. B 129, 365 (1977). 
[43] G. J. Alner et al., Z. Phys. C 33, 1 (1986) 
[44] K. Aamodt et al. (ALICE), arXiv: 1004:3034 (2010) 
[45] F. Abe et al., Phys. Rev. D 41, 2330 (1990). 
[46] W.D.Morse et al, Phys. Rev. D 15, 66 (1977) 
[47] V. Khachatryan et al. (CMS), JHEP02, 41 (2010) 
[48] V. Khachatryan et al. (CMS), Phys. Rev. Lett. 105, 

022002 (2010) 

[49] J. Klay et al, Phys. Rev. C 68, 054905 (2003); L. Ahlc 
et al., Phys. Rev. C 57, R466 (1998); B. B. Back et al., 
Phys. Rev. C 66, 054901 (2002); J. Dunlop, Ph.D. thesis, 
Massachusetts Institute of Technology (1999); L. Ahlc 
et al., Phys, Lett. B 490, 53 (2000) - the width a = 
•\/ln(s/4m 2 )/2 of the dN/dy distribution was estimated 
from Landau hydrodynamics, see e.g. P. Carrouthers and 
Ming Duong-van, Phys. Rev. D 8, 859 (1973) 
Phys. Lett. B 523, 79 (2001) 

[50] B.B.Back, Proceedings of Intl. Wkshp. Part. Mult. Rel. 
HI Coll. 1719 June 2004, Bari, Italy, J.Phys Conf. Ser. 5, 
1 (2005) 

[51] H. Heiselberg, Phys. Rep. 351, 161 (2001) 
[52] J. Benecke et al, Phys. Rev. 188, 2159 (1969) 
[53] B.B. Back et al. (PHOBOS), Phys. Rev. Lett. 94, 122303 
(2005) 

[54] I. Otterlund et al, Nucl. Phys. B142, 445 (1978) 



